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Abstract. The study of superfluid and Berezinskii-Kosterlitz-Thouless phases in 
exciton-polaritons requires an understanding of vortex dynamics in a dissipative 
unconfined condensate. In this article we study the motion of dynamic vortex- 
antivortex pairs and show that vortex pair stability defined as ordered motion 
as opposed to rapid separation or recombination is the result of balance between 
dissipative velocities in the condensate and interaction with thermal polaritons. 
The addition of a trapping potential is further shown to considerably enhance the 
lifetime of a single vortex pair in this system. These investigations have important 
consequences for interpretation of recent results and future investigations of two- 
dimensional superfluid phases in polariton condensates. 
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1. Introduction: Two-dimensional polariton condensation 

Superfluid behavior is well known in liquid helium and atomic Bose-Einstein condensates 
and extensive studies have been made of the properties of quantized vortices in these 
systems P E]. A newer type of condensate receiving a lot of current attention, 
occurs not in a gas or liquid of 'real' particles, but instead occurs in quasi-particles 
in semiconductors known as exciton-polaritons (hereafter referred to as polaritons) 
[31 S]. While the condensation of bosonic excitonic particles has been studied for 
many years, the strong coupling of excitons and photons reduces the effective mass 
by such a degree that condensation occurs in GaAs at temperatures of order lOK 
[3] and is expected at room temperature in other materials [5j. However, the small 
effective mass comes at the expense of finite particle lifetime, implying that although 
polaritons can reach equilibrium with the lattice [6j, a steady- state condensate is 
formed in place of a condensate in true thermal equilibrium. One of the advantages 
of a condensate embedded within a semiconductor matrix is the ability to manipulate 
the internal properties and dimensionality through material changes and the use of 
heterostructures. Though not strictly a necessary geometry, many of the recent 
observations of polariton condensation occur in two-dimensional (2D) quantum wells 
embedded within microcavities in order to reach the strong coupling regime. 

Superfluidity is an expected property of this condensate-like system [71 [8], and 
though not direct evidence of this, quantized vortices have already been observed pfTO]. 
Quantized vortices are topological excitations of phase coherent systems and exist 
commonly in response to an applied rotational field, but in two- dimensions spontaneous 
formation as a result of phase fluctuations is also possible [llj. Spontaneous vortex 
excitation is known to destroy long range order in a reduced dimensional infinite system 
preventing superfluidity |12[ at non-zero temperature. However, below a certain 
temperature, the thermal energy is not sufficient to generate distinct single vortices, but 
only the lower energy vortex- antivortex pairs (or vortex pairs for short). As the vortex- 
induced phase gradient is now largely localized between the vortex pair, superfluidity 
can be recovered. This 2D superfluid phase characterized by the presence of bound 
vortex-antivortex pairs is known as the Berezinskii-Kosterlitz-Thouless (BKT) phase, 
occurring at temperatures below Tbkt [14j. 

A Bose-Einstein condensate like transition can also be recovered in a 2D system 
providing the condensate is confined to a finite area of a size comparable to the phonon de 
Broglie wavelength in the fluid, thereby excluding long- wavelength fluctuations through 
a discrete density of states [15j. Although the transition from a normal (thermal) 
state to a BKT phase is discontinuous at Tbkt^ the transition from BEG to BKT is 
a continuous function of condensate size and temperature, both contributing to the 
vortex pair density in the system fT6] . The polariton condensate is thus an ideal system 
in which to study crossovers between various 2D condensate and superfluid phases. 

The BKT phase has previously been reported in superfluid ^He films [17j, atomic 
hydrogen films [I8] and dilute gas atomic BEG [I9]. In each of these cases however. 
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although convincing evidence for the BKT phase is presented, the microscopic nature of 
the state (the vortex-antivortex pairs) has yet to be observed. The smaU characteristic 
size scale of vortices (healing length ^) to the condensate size, results in a population 
of vortex pairs and/or free vortices far in excess of one, making single vortex pair 
observation difficult. The evidence supporting the BKT phase in atomic EEC however, 
is the observation of a proliferation of free vortices with increased temperature believed 
to be due to pair breaking when the temperature exceeds the pair binding energy. 

Recently however, direct evidence of a single vortex-antivortex pair has been 
reported in an exciton-polariton condensate [10] . The observation is based mainly on 
phase dislocations in interferometry experiments indicative of a single vortex-antivortex 
pair. Due to the low polariton effective mass m ^ 10~^m^ (m^ being the free electron 
mass), the vortices are significantly larger and for a typical condensate density and 
radius (L ^ 10/xm — 15/xm) one vortex pair can be comparable to the system size. 
Furthermore, unlike previous quantized vortex observations in polariton condensates, 
due to the low disorder in this sample {Vd ^ O.lmeV)^ this vortex pair is believed to 
be unpinned. This recent observation in particular suggests the necessity of a proper 
understanding of vortex and vortex pair dynamics in this system. 

In contrast to the superfiuid He and atomic BEC systems which generally have 
a particle number conservation over a dynamical timescale, the polariton condensate 
occurs in quasi-equilibrium, the result of a steady-state process of continuous stimulated 
scattering from a polariton reservoir (pump) and the finite polariton lifetime Tpoi ^ ps 
^20j . While the properties and dynamics of quantized vortices [211 122] and vortex pairs 
[23j have received considerable theoretical attention in superfiuid helium and atomic 
condensates, given its steady-state dissipative nature, the polariton condensate phase 
is distinct and the dynamical nature of quantized vortices presently unknown. The 
theoretical investigation of vortex dynamics presented here is directly relevant to the 
recent experimental observation of a single vortex pair and to further investigations of 
vortex nucleation and BEC-BKT crossover in a polariton condensate. 

2. Dissipative Gross-Pitaevskii equation 

Condensate dynamics are usually modeled using a form of the time dependent Gross- 
Pitaevskii equation (GPE) to describe evolution of the condensate order parameter 
'0(r, t) [24j. This technique has been extended to describe non-condensate components 
including quantum and thermal depletion and finite temperature effects [25l [26j . In this 
work we apply a dissipative GPE previously shown to contain the essential parameters 
necessary for simulation of polariton condensates [20l [271 E] • This dissipative GPE is 
coupled to a thermal reservoir population ni^(r, t) by stimulated scattering (i?(ni^(r, t))) 
and interactions (with coupling constant gR). Our study here specifically differs from 
other studies of vortex pair dynamics in that the non-condensate reservoir population 
is not uniquely determined by coupling with the condensate, but also by the spatial 
pumping profile. In addition because of the steady-state nature, this unconfined 
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condensate possesses a unique velocity profile. 

We start by defining the usual single particle Hamiltonian, 

Ho = -^ + VUr) (1) 

where m is the polariton effective mass, and Vext{^) is any external confining or disorder 
potential profile. The dissipative time-dependent Gross-Pitaevskii equation is then given 
as 

_ d^p{r,t) _ f ^ ih^^ ^^^^ ^ ^ 1^^^^^ ^ 2 



where the loss and gain terms of —ihjc/'^ and ihR{nR{r^ ^))/2 respectively describe the 
process of polariton decay as photons leaving the cavity at rate 7c (or alternatively 
with lifetime tc) and stimulated scattering into the condensate at a rate i?(ni^(r,t)) 
determined by the reservoir population distribution ni^(r,t). This reservoir population 
is described by a rate equation model, 

^^^^^ = Piiv)-^nnRir,t)-Rinn{v,t)mv,t)\' (3) 

where Pi{r) is the laser excitation profile (in photons//xm^]95) and is the reservoir 
polariton loss rate. 

Previous investigation of vortices with this model have not considered vortex pair 
dynamics, and it is not presently clear how a vortex pair should behave in a dissipative 
polariton condensate. Previous experimental observations of quantized vortices in 
polariton condensates have been carried out in materials with a level of disorder potential 
suflBcient to spatially pin the vortex, prohibiting any vortex dynamics and allowing 
simple experimental observation in time-integrated measurements. However, the recent 
clear observation of a spontaneously formed vortex-antivortex pair occurs in the GaAs- 
based system where the disorder is low enough for the vortex pair to remain unpinned 
and as such, dynamics and stability are critical to understanding this experimental result 
and is the main motivation for this research. 



3. Dissipative vortex pair dynamics 

The key parameters in the model of equations [2] and [3] which are expected to infiuence 
vortex dynamics in the absence of any trapping potential are the condensate density 
nc, reservoir density ur^ scattering rate i?(ni^(r,t)) and condensate polariton lifetime 
Tc. It is difficult to study the effect of any one of these parameter independently, as in 
this dynamic model changing either the lifetime or scattering rate both alter the ffnal 
steady state condensate population. As the condensate density critically determines the 
vortex size through the healing length parameter ^ = h/ ^J2mgcnc^> both R{nR{r^t)) 
and Tc are adjusted together in order to maintain an approximately constant steady 
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state condensate density. Equation [H describes the steady state relationship between 
these parameters (for a uniform system). 

ni^(r) Pi he- nc 

The scattering rate R{nR{r)) is given a hnear dependence on ni^(r) with coefficient 
Rsc The study of dissipative vortex- antivortex pair dynamics here specificaUy addresses 
Ga As-based polariton condensates, however general arguments will be drawn relevant 
to all other material systems. From a manually ascribed initial state, the condensate is 
allowed to evolve in time until a steady state distribution is achieved (constant energy 
and static particle distribution). The vortex pair phase profiles are then artificially 
imprinted along the y-axis (at positions ±dy/2 where dy is the vortex pair separation), 
with a mirror symmetry about the x-axis. The dynamics of an unpinned vortex- 
antivortex pair in a spatially-infinite conservative (non-dissipative) superfiuid is simple 
and well understood. Due to the interaction of the mutual phase gradients the vortex- 
ant ivortex pair undergoes a linear motion perpendicular to its dipole axis, with direction 
determined by the dipole orientation [2] (here, along the x-axis). The vortex pair will 
then travel at a velocity dependent on the vortex separation dy given by equation El 

where = l{h/m)z is the quantized circulation of each vortex with vector perpendicular 
to the x-y plane and / = ±1. In the absence of any dissipation, the vortex pair 
will maintain separation of dy. The introduction of a confining potential Vext{^) 
creates a boundary with which the vortex pair interacts, perturbing its one-dimensional 
motion. The form of the perturbation depends on the type of confining potential. In 
a conservative system, the vortex pair will track the boundary to preserve the energy 
of the vortex pair periodically returning to its starting position. Figure [TK shows the 
numerical solution of the conservative GPE in a square well trap (abrupt walls) giving 
an example of this type of motion. The trajectory in a harmonic trap is qualitatively 
similar, but this example is closer to that of a top-hat pumping profile considered in the 
rest of this paper for the dissipative condensate. 



3.1. Classification of vortex pair trajectories 

The vortex pair dynamics in a dissipative system are found to differ considerably. The 
initially imprinted pair with separation dy is found instead to choose from a continuum of 
different trajectories based on the parameters chosen for the simulation, some examples 
of which are shown in figure [lb- These motions range from the vortex pair splitting and 
leaving the condensate directly to recombining rapidly long before reaching the edge 
of the condensate. In these simulations, a condensate density nc ^ bYOfim''^ and a 
pumping spot of radius L = 15/xm are used. This condensate size is chosen such that 
the vortex pair is not initially perturbed by the presence of a boundary and the motion 
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Figure 1. The vortex pair trajectory for (a) a conservative condensate in a square 
well (abrupt walls of ImeV depth) circularly symmetric trapping potential of radius 
r = 15/im, and (b) dissipative condensates at various points in the parameter space 
depicting a range of paths between rapid separation and rapid recombination. The 
corresponding scattering rate coefficient is indicated (in iim?ps~^). The initial vortex 
pair positions are = ±4.5/im along the ?/-axis and the pumping spot radius is 
L = 15/im. 

largely independent of the condensate size can be initially observed, from which general 
conclusions are drawn. 

To gain an understanding of the effect of these parameters, a map of tir (size of 
circles) with variation in i?(ni^(r, t))/ni^(r, t) and tc is plotted figureO The background 
map (saturated at small scattering rate) is the steady-state analytical contribution 
(equation Hj) . In this parameter space the vortex pair trajectories are divided into two 
groups according to whether the vortex pair recombines in the center of the condensate 
before reaching the edge (depicted as blue circles) and vortex pairs splitting and leaving 
the condensate (red triangles). A clear transition between these two different groupings 
is apparent and observed to be largely dependent on i?(ni^(r, t)). 

If the initial vortex pair spacing dy is altered in this simulation, it does not affect 
this trend other than to slightly shift the transition between the two classifications of 
vortex pair evolution. This remains the case provided the vortex pairs are not initially 
overlapping slightly {dy < 2^) or the vortex is initially separated from the boundary by 
an amount of order \L — ry\ < 2^. In long timescale experimental scenarios where a 
statistical ensemble of initial vortex pair separations is expected this then corresponds to 
the same trend but with a blurred cross-over transition between the two classifications. 
The condensate population is also found to not impose much effect on this crossover 
(within experimentally reasonable ranges) as indicated in figure [2] where the different 
lines indicate different condensate population density. 

Attention is now turned to a map of vortex velocities constructed from this 
numerical data and shown in figure [3l The vortex velocities are presented as a fraction 
of the condensate sound velocity defined as c = \/gc^c/^- The distinction between 
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Figure 2. Map of ur (depicted as the circle/triangle marker size) with scattering 
rate R{nR{r)) and condensate polariton hfetime rc- The condensate density is held 
relatively constant in these simulations and corresponds to the different lines labeled 
with nc, but does not have a significant effect on the vortex trajectory. The vortices 
that leave the condensate are indicated as red triangles and those that recombine 
as blue circles. The background in this plot comprises the approximate steady-state 
analytical form of ur, with magnitude indicated by the colour bar in iim~'^. 



the two classifications of vortex dynamics is further clear in this data. For very low 
scattering rates (which corresponds to high reservoir population density), the vortex 
pair recombines rapidly. This vortex velocity drops rapidly from the sound velocity 
as the scattering rate is increased, stopping when the cross-over to radially-separating 
vortex pairs is achieved, remaining roughly constant thereafter at a fraction of the 
sound velocity. This analysis suggests two different mechanisms altering the vortex pair 
velocity vector. 
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Figure 3. A map of the numerical vortex velocities as a function of the R{nR{r)) and 
rc parameter space. The recombining vortex pairs are again indicated as blue circles, 
and the radially separating ones as red triangles. 
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3.2. Contributions to vortex velocity 

In a conservative condensate, the absence of any forces implies the vortex velocity 
will coincide with that of the local superfluid flow v^, = in correspondence with the 
familiar Magnus force (equation [6j) [281 [2]. rotating trapped condensate for example, 
the vortex lattice will rotate at the same angular velocity as the condensate and thus 
appear stationary in the rotating frame. In a dissipative system, i.e. through interaction 
with thermal population and excitations at flnite temperature, energy is transferred 
between these non-condensate populations and energetically unfavourable vortex states 
will decay by approaching the condensate boundary decaying via excitations at the edge 

fM = nci^ X (vl - Ys) (6) 

The superfluid velocity is deflned as being the superfluid velocity far from the 
vortex line. Thus, in the present vortex pair system, the superfluid velocity consists 
of velocity contributions from each of the two vortices (v^i and Vl2) and the radially 
dissipating polaritons due to repulsive interactions and a lack of conflnement. Neglecting 
the condensate boundary, to calculate the local velocity fleld at the vortex vl (at r^i 
with in the presence of vortex v2 (at r^2 with tvy2) and the radial velocity gradient 
(vc(r)) we can use, 

^-(^-i) = — r:^ — ;rT2 — + 

The second source of forces on the vortex pair is the interaction of vortices with 
non-condensate population which creates a drag force dependent on the interaction 
energy with this non-condensate population tir. This force is commonly broken up into 
longitudinal and transverse components. 



n BriR 



X (v^ - V^) 



where is the velocity of the normal component. The origin of these drag forces is 
usually attributed to interaction with thermally excited non-condensate modes, where 
the perpendicular component is commonly known as the lordanskii force [291 [30] . The 
form of the friction coeflicients B and B' have been evaluated theoretically |28l [3T] and 
experimentally [2j , though exact determination is very dependent on the exact details of 
the system and the contributions to the magnitude of the lordanskii force are particularly 
complex. We flnd however that conclusions can be drawn relatively independent of the 
exact magnitude of these coeflicients. 



3.3. Condensate velocity profiles 

As the condensate is unconflned, we expect there to be some radially dissipative 
condensate velocity in the system amounting to a continual radial loss of particles 
due to particle repulsive interactions, giving the form of vc(r). Figure shows 
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the numerically-evaluated steady-state unconfined condensate phase profile S{r). The 
condensate velocity is simply the gradient of the phase (v^ = h/mV S{y)) and the 
radial velocity profile can be calculated as in figure Hb. The point at which the 
velocity increases discontinuously is the approximate condensate boundary. Note that 
the velocities within the condensate possess a similar fraction of c as do those determined 
for splitting pairs in figure [3l 




Figure 4. (a) Steady-state unconfined condensate phase profile S{v) and (b) the 
radial cross-section of the velocity profile. The particle velocity within the condensate 
increases with the scattering coefficient R{nR{r)) /nR{r) in response to a reduction in 
drag forces through decreasing reservoir density nR{r) 



In the absence of energy dissipation, we expect vortices in the condensate to move 
with the same velocity vector as the local superfiuid fiow. The contributions to the 
magnitude of the vector vc{t) are from the same contributions as for the condensate 
sound velocity, namely, the condensate interaction energy gc\ip{T)\^ and the eflFective 
mass m. The condensate radial velocity also exhibits an inverse correspondence with 
the reservoir density nR{r) (induced by a reduction in scattering coefficient Rgc)- The 
presence of this thermal reservoir with which the condensate interacts introduces a drag 
force on the vortex pair. As the energy of a vortex pair is proportional to its separation 
(equation [9j) , in the presence of finite gRfiR^ energy will be transferred from the vortex 
pair to the thermal reservoir, inducing a drag force on each vortex towards the pair 
midpoint. 



Based on this simple analysis, the occurrence of two distinct modes of vortex motion 
differing from that of a conservative condensate appear to arise from contributions of 
two separate sources, namely the dissipative velocity of an unconfined condensate and 
the drag induced by interaction with thermal particles. Allowing the two forces to cancel 
(f^ + fi) = 0), the velocity vector vli of vortex vl is given by. 



BnR\K.yi\ , _ B'nR\K.yi\ 
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where kyi = tvvi/\f^vi\- Making substitutions for relevant for this specific case (see 
equation [71) reveals the contributions of these two sources to v^^i. If we assume the two 
vortices are initially at positions r^i = {dy/2)y and r^2 = —{dy/2)y with circulation 
vectors of tvyi = +{h/m)z and = —{h/m)z respectively yields the form of Vg for vl^ 

h 

v.(r,i) = + \Yc{d,/2)\y (11) 

Figure [5] shows a schematic of various contributions to the vortex velocity with 
(a) showing vc(r) and (b) and (c) showing the radially diverging and recombining 
vortex pairs respectively. The first term in equation [10] corresponds to the contribution 
of dissipative superfiuid fiow and the second term describing contributions of non- 
condensate interaction effects. While the exact trajectory depends on the accurate 
coeflBcients, it is clear that the first term contributes an outward velocity in the -\-y- 
direction with some curvature due to the +x component as a result of the vortex pair 
linear trajectory. Thus, on its own (i.e. when ur is small compared to n^), it describes 
a splitting and radially dissipating vortex pair as depicted in figure [5b. If we ignore 
the presence of the radial dissipative velocity {vc{y) = 0), the effect of the drag forces 
on vl in the second term of equation [10] can be described, which are proportional to 
nn. The perpendicular component labeled v^^^ = {BnR/2{nc + tir))^ x (v^ — v^) 
contains components in directions +x and —y while the longitudinal drag v^^n = 
— {B'nR/2{nc + nR))k x x (v^ — v^)] contains vector components in directions —x and 
—y as illustrated in [5)3. Thus, independent of the magnitude of the coefficients B and 
B' this set of equations generally describes a trajectory directed inwards (—^-direction) 
towards the mid point of the vortex pair. Clearly, if this analysis is performed again 
with the second vortex v2 the same results will be achieved with oppositely directed 
^-vectors, demonstrating that this simple analytical model agrees with our previous 
numerical results. 



(a) 1 


/ vc(r) \ ^ 


(b) 

i Tf 1 
o / 


(c) 

/ v^r& \ 

\ o ^ / 


•* — i — ^^^^^^^^^^^ 





Figure 5. Schematic illustration of the different components contributing to the vortex 
velocity, where the large dotted circle indicates the approximate condensate boundary 
and the small circles represent a sample of instantaneous vortex positions with time 
increasing to the right of the schematic, (a) shows the radial dissipative velocity vc(r), 
the result of repulsive interactions and lack of confinement. When the velocity in (a) 
is dominant, the vortex vector in (b) results and the vortex pair separates. When drag 
forces dominate due to large ur^ the vortex pair recombines as in (c). 
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3.4' Effect of trapping potential and interactions on vortex pair motion 

The application of a trapping potential will prevent radial particle escape by limiting any 
radial superflow present in the system and its presence can be used to check the effect of 
this superflow on vortex pair motion. Furthermore, the onset of the BKT-phase occurs 
not only with temperature, but also with conflnement size. In most current experiments 
the flnite laser pumping area is suflicient to restrict the condensate area (due to flnite 
lifetime and diffusion length). Experimental investigations of BKT transitions in a 
polariton condensate are likely to beneflt from the use of a trapping potential to conflne 
the gas, particularly if the lifetime is increased. Effective square well and harmonic 
proflle traps for polaritons have been previously experimentally demonstrated [32l [71 [33] 
allowing greater control over the condensate density proflle. 

In flgure[6]we show the numerical trajectories for (a) a diverging vortex pair and 
(b) for a recombining vortex pair. In flgure[6K, when a trapping potential is added, the 
vortex pair does not immediately split and leave the condensate, but recovers some of its 
circular motion remnant of the conservative scenario. This implies that the velocity of 
the radial dissipative flow as expected is the main contributor to perturbing the motion 
for this parameter range. In this case, the presence of a trapping potential extends 
the vortex pair lifetime by many times. In flgure [6b where the parameters dictate a 
rapidly recombining vortex pair, adding a trapping potential has no obvious effect on 
the trajectory, implying influence of any the radial condensate velocity is negligible for 
these parameters. 

Figure [6] also demonstrates the effect of turning off the interactions with the 
reservoir polaritons (setting gn = 0). In flgure [6^ the effect of turning off interactions 
only is to reduce drag forces towards the condensate centre and allows the vortex 
pair to leave the condensate more directly, via a shorter path. However, in flgure 
[6b for the recombining pair, turning off reservoir interactions completely prevents the 
recombination of the vortex pair and the vortices leave the condensate directly with the 
velocity vector comparable to Vg. 

4. Conclusion 

We have studied the dynamics and stability of a single vortex pair in a dissipative 
model of a polariton condensate. While vortex pairs are essentially stable against 
recombination and radial dissipation in conservative condensates, in an unconflned 
dissipative condensate, the vortex pair either splits and leaves the condensate or 
recombines quickly a short distance from the nucleation location. It is found that 
the cross-over of these two behaviours is a result of competition between a radially 
outward force due to the radially dissipating condensate polaritons and the interaction 
of the vortices with non-condensate population which strongly inhibits vortex motion 
and induces a drag force towards the pair midpoint as the vortex pair loses energy 
to the reservoir. We note that the long time-scale vortex pair stability appropriate 



12 



(a) splitting pair 



(b) recombining pair 



10 



3 

?3J 



-10 



▼ ▼ ▼ 9r = 
• • • trapped 
» « K untrapped 



10 



-10 



I • • trapped 
. M K untrapped 



tsIM 



-2 2 4 6 8 10 12 
X (/im) 



-2 2 4 6 
X {fim) 



10 12 



Figure 6. Vortex pair trajectories for parameters yielding (a) radially diverging vortex 
pair and (b) recombining vortex pair. Within these plots the effects of applying a 
trapping potential and turning off the interaction with the reservoir component are 
also displayed as separate trajectories. 



to conservative (atomic) condensates can be recovered in this system through the 
application of a trapping potential and pumping the system such that the reservoir 
density is not excessive. These observations thus have direct relevance to the 
interpretation of recent observations of a dynamic vortex-antivortex pair in a polariton 
condensate and to the extension of these studies to push the system into the BKT 
regime. 
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